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ABSTRACT 

Optimum  control  theory  is  applied  to  develop  a  guidance  law  for  homing 
missiles  for  which  the  airframe  dynamics  is  not  neglected.  This  work  extends 
the  applicability  of  the  laws  derived  in  a  previous  report  to  systems  character¬ 
ized  by  two  time  constants. 

The  effectiveness  of  the  optimum  controller  thus  derived  is  evaluated  in 
an  example  by  comparing  it  with  classical  proportional  navigation.  Additionally, 
extensive  quantitative  plots  of  the  controller  parameters  are  provided  as  an 
aid  in  specific  designs. 
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Section  i.  iNTRuDUCf  (0;; 


The  technique  of  proportional  navigation  (PN)  has  been  found  to  be  the 
most  satisfactory  method  of  guiding  homing  missiles  —  a  fact  established  by 
engineers  through  many  years  of  design  experience  rather  than  by  analytical 
proof,  in  PN,  an  attempt  is  made  to  mechanize  the  following  equation: 

y=NA;  2  <  N  £  5  , 

where  y  is  the  missile  angular  turning  rate  and  Ms  the  line— of-.^iu  angular 
rate  referred  to  inertial  space.  N  is  denoted  as  the  "navigation  ratio,"  and 
again,  the  range  of  acceptable  values  has  been  developed  mainly  through 
experience.  It  is  interesting  to  note  that  recently  it  has  been  rigorously  deter¬ 
mined  by  means  of  modern  control  theory  that  PN  is  indeed  optimal  in  that  fer 
unconstrained  control  effort,  the  miss  distance  at  intercept  is  minimized  in  the 
mean-squared  sense.  The  correspondence  between  PN  and  optimum  control 
(OC)  has  been  demonstrated  by  Bryson,  No,  and  Baron  [  1],  Janus  [2],  and 
Speyer  ( 3 J ,  among  others. 

A  previous  report  |4)  verified  the  optimality  of  PN  by  casting  the 
problem  as  a  "minimum  error  regulator"  problem  and  using  the  Ogata  1 5 ) 
method  of  solution.  The  weakness  of  the  cited  work  in  OC,  as  well  as  the 
classical  derivation  of  I’N,  is  that  the  plant  (missile)  is  assumed  to  respond 
instantaneously  to  guidance  commands,  i.c. ,  has  no  time  lags.  The  previous 
report  1 4 J  extended  prior  work  by  deriving  the  OC  for  a  missile  having  a  single 
time  constant.  It  was  shown  that  the  addition  of  the  time  lag  profoundly  affects 
the  steering  law  in  that  the  OC  lor  this  ease  requires  time  varying  gains. 

These  results  point  out  that  missile  dynamics  should  not  be  neglected  if  realistic 
OC's  arc  to  be  obtained. 

Typical  homing  missiles  cannot  realistically  be  characterized  by  a  single 
time  constant;  two  or  even  three  lags  are  required  for  adequate  modeling.  The 
purpose  of  this  study  is  to  derive  the  OC  for  a  two-time  constant  system,  and  to 
analyze  the  performance  of  this  controller  relative  to  conventional  PN.  The 
effect  of  using  such  a  controller  in  a  system  characterized  by  three  time  con¬ 
stants  will  also  bo  investigated;  it  is  conjectured  that,  although,  it  is  suboptimal 
for  this  case,  such  a  scheme  would  nevertheless  be  superior  lo  PN,  since  the 
latter  takes  no  time  lags  into  account.  The  method  of  solution  used  is  to  cast 
the  problem  in  a  form  such  that  the  Ogata  method,  which  has  been  described  ami 
verified  f  4  J ,  can  be  used. 


1 


Section  II.  THE  MINIMUM  ERROR  REGULnTorv  rrtOBLEM 


Any  linear  dynamical  stem  of  the  nVn  01  ler  can  bo  expressed  either  as 

an  order  differential  equation  or  a  set  of  n  first-order  differential  equations. 
The  latter  is  known  as  the  state  formulation  and  is  used  herein,  since  it  lends 
itself  to  matrix-vector  notation  and  manipulation.  It  will  be  assumed  that  the 
system  differential  equation  is  given  by: 

x  j=  Ax  +  Bu;  x  (0)  =  C  , 

where 

x  =  n  dimensional  column  state  vector 
u  =  r  dimensional  control  vector 
A  =  n  x  n  matrix 
B  =  n  x  r  matrix 

and  where  the  following  index  is  to  be  minimized: 

T  T 

J(C,T)  =  x*(T)Px(T)  +  /  x*(t)Qx(t)dt+  /  u*(t)R(t)u(t)dt  . 

- - 0  - - • - -  0  - — - 

Terminal  State  Control  cost 

state  weighting 

weighting 

The  *  symbol  denotes  the  conjugiite  transpose  of  the  vector,  or  simply  the 
transpose  for  real  vectors,  and  P,  Q{t),  lift)  are  matrices  of  appropriate 
dimensions. 

Ogata  shows  that  the  optimum  controller  for  such  a  system  can  be 
obtained  by  solving  the  nonhomogeneous  matrix  Riceati  equation: 

-SA  -  A*S  -t-  SBR"‘  (0)B*S  =  Q(0)  . 

If  the  matrices  A  and  B  are  constant,  i.e. ,  if  the  system  is  stationary,  the 
above  matrix  equation  can  be  solved  in  closed  form  for  the  time  T: 

S(T>«  {[<PU(T)  +  <^>22  (T>P|  {  0n(T)  +  «0(T)P]-‘}  , 
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where  the  0  are  obtained  from  partitioning  the  matrix 


ajui  M  is  defined  as 


M  = 


A 

4>,t<  T> 

1 

l  ^>,2  (T ) 

<p2im 

!  022  <T) 

r 

L 

-A  1 

1 

BR_1(0U3* 

Q  (0)  ! 

A':' 

from  which  it  can  be  seen  that  M  is  known  from  the  problem  statement  and  the 
performance  index. 

Once  S(T)  is  known,  the  optimum  control  vector  can  be  obtained  from 
the  expression 


uopt(t)  =  -F(T-t)x(t,  . 

where 


F  (T  -  t)  =  R-1(t)B*S(T  -  t)  . 

In  block  diagram  form,  the  OC  can  be  depicted  as  in  Figure  1. 


FIGURE  1.  OPTIMAL  CONTROL  SYSTEM 


The  solution  of  the  subject  problem  can  then  be  summarized  as  follows: 

1)  From  system  state  equations,  and  the  given  performance  index,  the 
following  matrices  are  known: 

A,  B,  P,  Q,  R 

2)  The  matrix  M  can  be  formed  from  above 

I 

-A  I  BR-'(0)B* 

M  =  “""■}■ - 

q  (0)  ;  A* 

MT 

3 1  From  knowledge  of  M,  e  can  be  found.  There  arc  several 

methods  for  computing  this,  but  the  use  of  the  Laplace  transform  is 
often  the  most  convenient.  This  method  involves  the  relationship 

Jil[eMTJ  ~  [si  -  M  ]  _1  , 

where  s  is  the  Laplace  operator  and  I  the  unit  matrix. 

MT 

4)  Once  e  is  known,  all  of  the  <f>^  are  obtained  from  the  relation 

<Pu(T)  \  0,2  (T) 

* 

-MT)  !  022  (T) 

From  this,  S(T)  can  be  computed: 

S(T)  =  {[02i(T)  +  ^(T)P][0„(T)  +  M^P]-1}  . 

5)  Once  S  (T)  is  kn<  vn,  u  (t)  and  F  (T  -  t)  are  immediately 
obtainable: 

F(T  -  t)  =  R-  j(t) B^S  (T  -  t) 


MT 

e 


u  (t)  =  F(T  -  t)x(t)  . 
opt 

The  optimum  system  thus  mechanized  will  minimize  the  given  quadratic  index  in 
the  solution  interval  0  ^  t  ^  T.  It  should  be  noted  that  the  optimum  solution  is 
given  in  terms  of  tin;e-to-go  (T  -  t)  rather  than  elapsed  time  t. 
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Section  iii.  OP  l  iMUM  CONTROLLER  FOR  A  TWO*TIME  CONSTANT  MISSILE 


The  conventional  overall  loop  can  be  depicted  as  shown  in  Figure  2  and 
can  bo  readily  rearranged  with  target  acceleration  (  y  )  as  an  input  in  the 

manner  of  Figure  3.  By  assuming  two  time  constants  for  the  dynamic  lag  of  the 
missile-autopilot  combination,  the  model  of  Figure  4  is  obtained.  This  figure 
also  includes  an  exponential  decay  model  for  target  acceleration,  to  the  left  of 
the  dashed  line  in  Figure  4.  The  contribution  of  the  target  acceleration  has, 
however,  been  previously  computed  [4]  and,  since  it  remains  invariant  with 
the  order  of  the  plant,  it  does  not  need  to  be  explicitly  considered  herein.  With 
reference  to  Figure  4,  the  problem  can  be  staled  as  follows:  Given  the 
observable  states  x,,  x2,  x3,  and  x, ,  the  control  vector  is  determined  that  will 
minimize  the  miss  distance  at  intercept  only,  i.e. , 

y  (t)|  =  minimum  in  the  mean  squared  sense  subject  to  a  constraint 

1 1  on  available  control  effort. 

The  state  equations  can  be  written  directly  from  Figure  4;  the  target 
acceleration  term  is  ignored  for  now  and  will  be  appended  later,  since  ft  is  the 
same  as  for  the  cases  previously  treated. 

*1*  *2 
X2  =  -  g  X3 

X3  -  -  ^2  X3  +  w2 
X4=  -  wix4  +  w,n  . 

L 


I# 


FIGURE  2.  TYPICAL  HOMING  BLOCK  DIAGRAM 
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FIGURE  4.  AUGMENTED  HOMING  SYSTEM 


In  the  vector-matrix  form,  these  become 


0 

1 

0 

0 

*1 

0 

0 

0 

-  g 

0 

X, 

0 

0 

0 

-  cr. 

^2 

*3 

f 

0 

0 

0 

0 

-  u, 

x4 

“l 

The  generalized  index  J 

T 

J  =  x*(T)Px(T)  +  f  |x*(t)Q  (t)x(t)  +  u*(t)H  (t)u(t) ) dt 

'o 

is  reduced  to  the  desired  form  by  defining 

0  0  0 

0  0  0 

0  0  0 

0  0  0 


;  Q"|0)4x4;  u=*=scalar 


which  yields  the  form 


J  =  Xj(T)2  +  X  J  n  5dt  . 


ml  1  •  1  •  t  _  .  1  -4.  3..  X - - 4.  _  J  - - J  rT~I1 - *  £  -- 

x  ne  aix»ve  inuex  ini!!iiiiiz.c*&  univ  v  .  eu.  ititei  cebi  utantu,  me*  mmi  ia  mi  a» 

d 

obtained  next: 


from  v.hich  [si  -  M]  can  be  formed; 


Since  the  above  matrix  is  almost  triangular,  it  is  most  easily  inverted 
by  reducing  the  augmented  matrix, 

[si  -  M  i  I]  , 

to  the  upper  triangular  form.  The  inverse  is  obtained  by  solving  the  triangular 
matr  ix  as  a  set  of  algebraic  equations  using  each  column  of  the  augmented 
section  in  turn.  This  process  yields  the  matrix  [si  -  M]-1  as  shown; 


8 


(85-Wj)  (8+Wj) 


s(»+-o),)  (8+Wj) 


(S+W2) 


Correctness  of  this  matrix  can  be  verified  by  multiplying  it  by  [si  -  M  J  and 
observing  that  the  unit  matrix  results.  The  0  needed  for  the  compilation  of 
S(t)  are  now  available,  since  J 


|sl  -  MJ 


■!  _ 

02,(S)  <?>22(S> 


022(s)  | 

The  computation  of  S(t)  requires  the  computation  of  the  following  two  matrices  - 

I  4> 21<t)  +  <^22  ft )  F  i 

<t)P]  . 

Since  P  is  a  number  matrix,  and  the  addition  operation  is  invariant  in  the  inverse 

Laplace  transformation,  the  terms  within  the  brackets  can  be  computed  prior 

to  determining  the  inverse  transform  of  the  d>...  Also,  not  all  of  the  terms  of 

ij 

each  0  will  be  needed;  therefore,  the  matrix  elements  will  be  denoted 

symbolically.  This  reduces  computational  complexity  considerably.  The 
matrix  <f>21(t)  is  the  null  matrix;  therefore  the  quantity  [  <A?)(t)  +  <f>22(t)P]  is 
merely  the  transform  of  the  first  column  of  <J>22(S),  to  be  denoted  symbolically' 
as 


|cf,21(t)  +  4>22(t)P]  = 


0  0 
0  0 
0  0 


0  0 
0  0 


0  0 


The  form  of  [  +  4>eP)  is 


<J>  1 1  +  <p  12  P 


an 

a  (2 

a  is 

a  14 

»2I 

a22 

a  23 

a24 

a31 

0 

a33 

a34 

a4  1 

0 

0 

a44 

Inversion  by7  partitioning  into  2x2  matrices  is  simplest  in  this  case;  given  a 
matrix 


A  j  B 

s  =  —  r - 

C  l  D 
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and  ita  inverse 


K 

M 


L 

N 


The  relationship  between  the  two  matrices  can  be  established  with  the  aid  of  the 
auxiliary  relations, 


[ s"~ isJ  =  [ss  tJ  =  I 

which  yield  the  equations 


K  =  (A  -  BD~  'C )  “ 1 

L  =  -KBD'1 

M=  -  D"'CK 

N  =  -  D-  1  -  D' 

The  application  of  the  above  formulas  to  the  symbolic  representation 
of  [<pu  +  <pi2 PJ  yields  the  desired  inverse; 


aMaSJa4l  “a|?a3Salj  'a22a13utl+al2a23a44  ~a  Ma22a33+a  13llJ4a22"a44a23a  U+  a  12a24aJ3 


where 


^*i=  fa ; ;Sz;?-.iia44  "  S2.a13‘M4“3!  +  a22al3*13.iil4  J 


a22H14a4)ajS  a12a21“33a41 


+  aEa23a44a3j  -  ?-t2&Ua23aU+  a  I2a24a33a4  1)  * 


It  will  be  sliovra  that  only  the  first  row  of  1  0n+  12 P ] ” 1  is  needed.  For 

this  reason,  the  1  jmaininf.  rows  were  not  computed.  The  matrix  S(tj  can  now 
be  computed  since  S (t)  =  [$2]  +  <p22 P)  f$n+  ^PJ”1; 


U 

it 

K  j 

h  Kj. 

h  K3 

h  K4 

i 

Ks 

i  K, 

i  k3 

i  k4 

j 

Ki 

j  k2 

j  k3 

j  K4 

k 

K, 

k  K2 

k  K3 

k  K4 

12 


■(kil 


i  wr>  %  i — k-»  mai'mmit'V- 


where 


K|=  a22a33aill 
K2  3  _ai2a33a44 

K3  3  _a:2aI3a44+  a12a2Sa44 

K4  =  ”a  14a22a33  +  a  13a34a22  ~  a34a23a12  +  a  I2a24a33  • 

Finally,  the  optimum  controller  can  now  be  written  from  the  previously 
stated  formula  F  (T-t )  =  R" 'B  :!S  (T-t ) 

F(T-t)=^  lK,K,  K,X,| 

go 


and 


u  ,  A  \  =  -  F/t  \  xft) 
opt  {  go  J  \  go ) 


The  remaining  task  consists  of  obtaining  the  quantities  k,  C,,  K K2,  K3, 
and  K4  as  a  function  of  the  specific  problem  variables;  for  the  subject  problem 
this  is  algebraically  an  onerous  task.  The  first  step  consists  of  finding  the  a 


and  k  by  inverse  Laplace  transforming  the  elements  of  (4>n+  <pn P) .  The 
details  of  this  computation  are  relegated  to  Appendix  A.  Knowledge  of  the 
values  of  k  and  the  a^.  allows  Cj  and  the  K^,  and  consequently,  F(t)  to  be 

computed.  By  omitting  the  algebraic  details,  (t  —  t  is  implied) 

t?o 


ij 


-  g 


-co,t 


>  r. 


co. 


W  e’"*4 
+  co2  (co,-co2) 


tui.  +  co„H 


IX  +  g2  w,  <1)1 


is  obtained,  where 


I  = 


2(co,t— 1)  —  co.t  2  (cOjt-1)  — co,t 

— s— 4 - -  e  -  — 5 — —  e  t 

a>j  (corco2)  co2  (tOj-Wj) 


-  CO, 


2co,s(corco2)5 


it)  ?  — 9fiir  f 

— 5 - L  a  G  ~  z  ~ 

((1). 2 


2 CO 2  (CO  j-  Ct^ )  ^ 


(C0,2-C022)  (C0,-C02) 


CO]2 


2co7r(a’i‘-°V!)  2C0,-’  (co^-Wj2) 


CO,2 
3  1  2  \ 


(COj  +  U>&)(u)tt-1)  (C02t~l)  t3 
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It  is  interesting  to  note  that  the  terms  Gj,  G2,  and  G3  are  identical  to  the 
three  equivalent  terms  obtained  for  the  one-time  constant  system  [4].  It  should 
also  be  noted  that  the  terms  of  G3  and  G4  within  the  parentheses  are  all  identical 
in  form.  This  is  significant  in  that  the  same  algorithm  can  be  used  to  compute 
these  terms. 

The  equations  for  N  and  G4  are  not  valid  for  ce,  =  since  the  term 
-  w2)  becomes  zero.  This  difficulty  is  resolved  by  computing  a  new  set  of 
equations  for  those  terms  valid  only  for  the  case  where  Vj-  By  appropriately 
combining  the  terms  of  N  and  Gt,  the  elements  containing  (w,  -  u>2)  can  be  iso¬ 
lated  into  the  indeterminate  form  0/ 0  when  a>(=  u>2.  LTIospital's  rule  can  then  be 
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applied,  where  the  differentiation  is  performed  with  respect  to  either  to.  or  u'2. 
With  the  indeicrniinaele.8  removed,  the  substitution 

a' ,  =  a.',  -  w 


is  made  into  the  resulting  equation.  Application  ot  this  technique  results  in 
the  equation 


N 


i  •)  /  -cat 

gJ  t2  ^tot  -  2  +  cut  e  +2  c 


-wt' 


toj) 


X  +  g’  to  1 

<«1  =  “V 


where 


(to,=to2) 


*  a  —  Cc’t 

r  _3_  2foe _ _  4te 

3  "  4  to*  ”  a' 


G  =  -  — 

4(Ci.'1-  a'2}  C*r 


or 


4c 

-+ — r 

CO 


—cot  a  -2o.'t  -2  cot 


tJe 


5te 


13e 


-2wt 


2to 


2  to 


2tf  4t 
Tar"-  ca  er 


’  / 

Wt  tot  \ 

-tot 

[ 

2  +  tot  e  -  2o 

+  -We 

[\ 

tot  ?  ) 

e  t  ' 

1  +  t 

We  now  have  equations  that  are  valid  for  all  ir  .  If  a>2  is  allowed  to 

approach  infinity,  i.e. ,  the  second  time  lag  is  eliminated,  the  above  equations 
should  and  do  reduce  to  the  ones  previously  derived  (  4)  for  the  single  lag  system. 


or 


The  optimum  controller  can  now  be  written 

r  x'~ 

I  x2 

u  -  n  —  N  |G.  G-/  Gj  C4]  I 
Opt  c  ‘  X3 

L  *4  . 


n  =  N  I G  j  x «  +  Gj  k2  +  Gs  x3  +  G4  x1]  . 
c 

However,  it  is  seen  (Figure  4)  that  Xj-  y^  and  x2  =  y  .  Therefore,  in  terms  of 
these  parameters,  and  the  derived  values  of  Gj  and  G2,  n  can  be  rewritten 


n 


=  N 


\6  TP 


+  N  G2  +  K  Gi  x4  . 


The  reason  for  this  representation  is  that  y  and  yf}  are  not  usually  explicitly 
measured,  since  conventional  seekers  measure  the  line-of-sight  rate  \j  since 
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A  =  y  /H,  where  R  is  the  missile-target  range,  the  following  relationship  holds: 

/. 


V  t 
c  go 


v  r" 

\  go 


where 


V  =  closing  velocity  =  -R. 

By  use  of  the  above,  the  control  equation  can  be  rewritten 


NV 

n  = - -  A  +  NG3  x3  +  NGj  x,  . 

eg  J  44 

The  first  term  on  right  side  of  above  equation  is  recognized  as  being  analogous 
to  classical  PN,  except  that  N  is  time-varyirg  instead  of  fixed.  The  optimum 
controller  also  has  two  additional  feedback  components  from  the  states  x3  and  x4. 

The  inclusion  of  the  target  acceleration  model  of  Figure  4  requires  one 
additional  term  in  the  control  equation,  the  latter  having  been  derived  [4].  The 
controller  would  take  the  form 


NV 

ne  ~  ^  +  NGs  *3  +  NGi  x 4  +  NGs  *6  ’ 


Again  it  is  noted  that  G<;  has  the  same  form  as  G3  and  the  components  of  G4. 

All  of  the  analytical  parameters  required  to  construct  the  optimum  con¬ 
troller  have  been  obtained.  In  the  next  section  a  quantitative  study  of  a  system 
using  such  a  controllei  il  be  performed. 
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Section  !V„  QUANTITATIVE  ANALYSIS  OF  A  TYPICAL  HOMING  SYSTEM 


I.  Analysis  of  the  Controller  Equations 

The  optimum  homing  system  requires  computation  of  N,  G3  and  G4. 
At  first  glance  this  appears  to  be  a  formidable  task  since  these  quantities  are 
elaborate  functions  of  t  ,  Wj,  and  w2.  However,  for  a  given  missile,  and  w2 

would  remain  relatively  constant,  at  least  for  a  limited  altitude  band;  in  this 
study  they  will  be  assumed  to  be  constant  and  the  required  functions  thus 
become  dependent  only  on  t  for  a  given  case.  The  more  general  case  of 

varying  u)j  and  w2  will  be  discussed  later. 


Due  to  their  complexity ,  it  is  in  most  instances  impractical  to  compute 
the  required  parameters  in  real  time  from  the  given  equations.  A  search  for 
simpler  mechanizations  requires  first  that  the  nature  of  the  subject  functions  be 
studied.  For  this  purpose,  values  of  N,  G3  and  G4  have  been  computed  for  a 
wide  range  of  values  of  w,  and  w2  by  using  a  digital  program.  Plots  of  these 
functions  are  shown  in  Appendix  B.  It  is  seen  that,  although  they  originate  from 
complex  equations,  the  plots  of  N,  G3  and  G4  show  these  quantities  to  be  w'ell 
behaved,  continuous  functions  oft  .  Thus,  for  a  given  configuration  where  u>, 

and  u2  are  constants,  these  quantities  could  be  easily  generated  as  functions  of 
t  by  analog  (diode  function  generator)  or  digital  (read-only  memory) 

techniques. 


2.  Comparative  Evaluation  of  Optimum  Controller  and  Proportional  Navigation 

The  increased  complexity  of  the  OC  over  PN  indicates  that  to  be 
feasible  the  OC  must  yield  a  considerable  payoff  over  PN  in  terms  of  reduced 
miss  distance.  For  this  reason,  a  comparative  evaluation  using  simulation 
techniques  was  performed  on  a  simple  model,  using  both  the  OC  and  PN 
concepts. 

The  sample  problem  concerns  a  small  homing  missile  in  a  low-altitude 
air-to-ground  engagement,  with  the  target  assumed  stationary.  The  pertinent 
geometry  is  depicted  in  Figure  5.  Ihe  model  of  Figure  4  with  the  OC  appended 
can  be  formally  represented  as  shown  in  Figure  6-a.  However,  since  conven¬ 
tional  seekers  measure  A,  a  more  convenient  representation  is  that  of  Figure 
6-b.  This  model  v/as  used  to  generate  two  digital  simulations  using  the  CSMP 
language.  The  block  diagrams  of  these  pregrams  are  shown  in  Figures  7-a 
and  7-b  respectively;  the  first  of  these  is  a  direct  simulation  of  the  system, 
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FIGURE  5.  SIMPLIFIED  HOMING  GEOMETRY 


and  the  second  is  an  adjoint  model  used  for  study  of  stochastic  c  ffects.  Each 
model  was  used  in  the  PN  mode  where  N  was  held  to  a  constant,  and  in  the  OC 
mode  with  time  varying  coefficients.  The  specific  parameters  are  summarized 
in  Table  I  for  each  case. 


TABLE  I.  SIMULATION  PARAMETERS 


Par;  neters 

PN 

OC 

1 

1 

uj2 

2 

2 

N 

4 

* 

C3 

0 

* 

0 

* 

R 

o 

10,  000  ft 

.0,  000  ft 

V 

1,  000  ft/ sec 

1 ,  000  ft/  sec 

*  Computed  as  a  function  of  t^  from  OC  equations. 

The  complexity  of  the  equations  for  N,  Gs  and  C4  precluded  "on  line" 

computation.  These  values  were  precomputed,  and  entered  in  the  simulation 

as  tables.  The  CSMP  function  generator  takes  these  tabular  values  and 

performs  linear  interpolation  for  values  of  t  lying  in  between  table  entries. 

go 

No  effort  was  made  to  accurately  simulate  the  subject  parameters  for  very  short 

ranges.  Due  to  mechanization  limitations,  most  seekers  have  a  "blind  range" 

of  several  hundred  feet,  and  are  unable  to  accurately  measure  A.  Therefore, 

N,  Ci  and  G a  were  linearly  decremented  to  zero,  between  t  =0.5  and  t  =0 
i  4  go  go 

sec.  These  approximate  functions  are  shown  in  Figures  8-a,  8-b,  and  8-c. 

The  simulation  study  directly  evaluated  PN  with  a  gain  of  4  to  the  approxi¬ 
mate  OC  using  the  parameter  approximation  described  above.  A  number  of 
input  perturbations  were  considered  and  the  resulting  miss  distance  was 
measured.  Specific  disturbances  were  as  follows: 

Case  I:  Impulsive  heading  disturbance  of  1  degree  occurring  at  an 
arbitrary  t 

go 
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FIGUiiE  8-c.  SIMULATED  NG4  FUNCTION 


Case  II:  Impulsive  disturbance  on  A  of  l/G  degree  per  second  occurring 

at  an  arbitrary  t 
'  go 

Case  III:  Additive  wide  band  noise  on  A  with  flat  power  spectrum  in 
region  Q  -  20  hertz 

Case  IV:  Additive  narrow  band  noise  on  A  with  a  single  peaked  power 
spectrum  as  shown  in  Figure  9. 

The  simulation  results  are  shown  in  Figures  10  through  13,  respectively, 
for  each  case.  It  is  evident  that  a  dramatic  reduction  j.n  miss  distance  results 
from  use  of  the  OC;  thus,  the  added  complexity  of  implementing  such  a  control 
scheme  may  in  certain  instances  be  well  worth  the  effort. 


I 


I 


The  OC  equations  are  optimum  only  if  the  plant  has  two  time  constants. 
The  effect  of  implementing  this  OC  in  a  higher  order  system  was  also  studied. 

A  third  time  constant  of  i  second  was  added  to  the  plant,  and  this  system  was 
again  compared  with  the  use  of  a  constant  N.  These  results  are  shown  in 
Figures  14  and  15,  Evidently  the  OC  is  still  superior  to  PN,  although  perform¬ 
ance  has  degraded  somewhat.  This  is  a  significant  result  in  that  most  practical 
systems  are  of  high  order.  The  above  indicates  that  if  an  optimum  scheme  is 
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derived  that  accounts  for  only  the  predominant  time  constants  (two  in  this 
report),  the  resulting  system  performance,  although  suboptimul,  would  still 
be  superior  to  PN. 


3.  The  Effective  Navigation  Ratio 

In  classical  PN  the  gain  N  denotes  the  static  ratio  of  missile  turning 
rate  to  line-of-sight  rate.  The  gain  N  is  usually  a  small  number  between  3  and 
5.  The  plots  of  N  in  Appendix  B  show  that  for  the  OC,  N  takes  on  very  large 
values  relative  to  the  PN  ease.  One  is  tempted  to  conclude  that  the  OC  requires 
much  higher  turn  rates  than  PN.  This  is  not  true,  however;  in  Figure  6  it 
can  be  seen  that  the  two  nonunity  feedback  terms  contribute  to  the  static  gain 
of  the  loop.  The  actual  value  of  y/A  can  be  computed  by  reducing  the  two  inner 
loops  and  obtaining  a  single  equivalent  expression  for  these.  The  resulting  gain, 
ignoring  missile  dynamics,  is  the  navigation  ratio  for  the  optimum  system. 

This  is  given  by 


N 


N 


off  1  +  N(G3+G4) 
large  dynami 

considered  herein  it  hovers  in  the  region 


N  does  not  have  a  large  dynamic  range.  For  the  numerical  example 
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TIME  BETWEEN  APPLICATION 
OF  DISTURBANCE  AND  INTERCEPT  (sec) 

FIGURE  10.  MISS  DUE  TO  RANDOM  HEADING  DISTURBANCES 


TIME  BETWEEN  APPLICATION 
OF  DISTURBANCE  AND  INTERCEPT  (sec) 

FIGURE  11.  MISS  DUE  TO  RANDOM  DISTURBANCE  ON  A 
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2  <  N  <  3; 


thus,  the  OC  does  not  require  increased  missile  performance. 


4.  The  Effect  of  the  Parameter  X 

The  quantity  A  represents  a  weight  on  the  parameter  N.  For  large 

values  of  t  it  1ms  little  or  no  effect,  but  reduces  N  heavily  when  t  is  small, 
go  •’go 

This  effect  can  be  seen  in  Appendix  B,  if  one  compares  the  plots  of  N  for  X  =  0 

with  those  for  A  =  1.  Since  X  reduces  N  near-  intercept,  it  obviously  reduces 

required  missile  maneuver,  since  N  ^  is  also  reduced  in  magnitude.  Since  X 

is  thus  a  significant  design  parameter,  how  does  one  choose  a  "good"  value? 
This  problem  is  intimately  tied  to  each  specific  missile  design,  since  it  relates 
directly  to  available  maneuver  capability.  From  Figure  6-a,  the  called- for 
maneuver  in  G  can  be  written 


n 

c 


N 

1  +  N(G3  +  G4) 


V  A 
c 

g 


eff  ■ 


By  defining  the  maximum  available  maneuver  as  u  ,  the  following  expression 
can  be  derived  from  the  above  equation: 


N  =  . 

eff  V  A 
c 

The  terms  on  the  right  side  of  this  equation  can  be  established  for  a  given 

missile  and  engagement  configuration.  Thus,  the  maximum  practical  value  of 

N  ,,  is  defined,  and  this  in  turn  infers  a  specific  A,  since  there  is  no  need  in 
eff 

generating  values  of  that  result  in  a  called-for  missile  maneuver  in  excess 

of  performance  limits.  The  penahy ,  of  course,  is  that  for  A  >  0,  the  miss 
distance  will  increase  over  that  obtainable  with  X  =  0.  This  is  illustrated  in 
Figure  16  for  the  example  system. 


29 


5 


FIGURE  16.  EFFECT  OF  VARYING  \  ON  MISS  DUE  TO  RANDOM 
DISTURBANCE  ON  A 
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Section  V.  MECHANIZATION  CONSIDERATIONS 


A  mathematically  elegant  solution  is  useless  if  hardware  implementation 
is  not  practical.  As  mentioned  before,  N,  G3  and  G4  can  be  generated  by  analog 
(diode  function  generators)  or  digital  (read-only  memories)  means.  The  input 
to  such  devices  would  be  t  ,  thus  requiring  a  measurement  of  this  quantity. 

For  stationary  or  slow-moving  targets,  t  can  be  estimated  closely  by  know¬ 
ledge  of  the  initial  range  and  missile  velocity  profile.  Intercept  of  fast  targets 
that  contribute  considerably  to  the  closing  velocity  is  more  difficult,  since 
ranging  must  be  obtained  during  flight.  This  requires  tracking  of  both  inter¬ 
ceptor  and  target  by  either  ground  radar  or  an  active  ranging  device  on  board 
the  interceptor. 

Another  major  consideration  is  that  the  OC  requires  that  all  system 

states  be  measured.  The  actual  missile  acceleration  can  be  measured  with  an 

accelerometer;  y  and  y  are  implicitly  measured  by  the  seeker  as  previously 
d  d 

shown.  The  remaining  state  (x4  in  Figure  4)  cannot  be  measured  if  the  two 
time  constants  are  assumed  to  represent  the  airframe.  However,  if  in  an 

n^1  order  system  (n-1)  states  are  measurable,  the  one  can  be  estimated  by 
"observer  pole"  techniques  |6~8] .  Finally,  if  the  target  acceleration  term  G5 
is  to  be  included,  y  must  also  be  measured.  This  requirement,  although  diffi¬ 
cult  to  meet,  is  no  different  from  the  requirement  for  implementation  of  biased 
PN,  where  a  bias  term  is  appended  in  order  to  reduce  the  miss  caused  by  target 
evasion. 

Certain  types  of  missiles  are  required  to  operate  through  a  wide  altitude 
band,  where  air  density  considerably  affects  the  airframe  time  constants;  for 
such  missiles,  Wj  and  <j2  cannot  be  assumed  constant.  This  problem  also  exists 
in  conventional  systems  required  to  operate  in  such  an  environment.  The 
solution  is  to  generate  a  family  of  functions  N,  G3  and  G4,  each  valid  for  a 
specific  altitude  band.  The  controller  would  then  be  "band-switched"  to  the 
appropriate  functions,  dependent  on  the  altitude. 

The  values  of  N  shown  in  Appendix  13  are  plotted  in  log-log  coordinate 
paper.  It  is  seen  that  for  the  vast  majority  of  eases,  N  can  be  approximated  by 
a  few  straight  line  segments  in  the  log-log  plane.  This  suggests  a  very  simple 
mechanization  scheme  for  N.  The  above  linear  relationship  implies  that 

in  N  =  C.  +  b,  ir.  t  , 

1  1  go  .  . 

region  1 


where  is  a  sure  able  intercept  and  is  the  signed  slope  of  the  straight  line 
approximation. 

A  family  of  equations  as  above  can  be  mechanized  as  shown  in  Figure  17. 


Section  VE.  CONCLUSIONS 


The  problem  solved  herein  is  an  extension  of  the  work  begun  in  1968  ( 4  ] . 
The  basic  homing  system  has  been  augmented  by  two  time  constants  and  cast  as 
a  minimum  regulator  problem.  Simulation  using  the  resulting  controller 
equations  indicates  that  the  optimal  approach  results  in  a  guidance  scheme 
superior  to  PN. 

Practical  problems  have  been  considered  and  it  is  concluded  that  the  OC 
concept  is  feasible  from  a  mechanisation  viewpoint.  The  added  complexity 
must  be  weighed  against  the  substantial  payoff  in  drastically  reduced  miss 
distances. 

The  behavior  of  the  OC  derived  herein  is  consistent  with  previously 
derived  results.  For  very  large  t  ,  N  is  assymptotic  to  3,  As  the  time 

constants  approach  0  (w  — ■  «)  N  approaches  3  and  Cl3  and  G4  approach  0.  Thus 

for  very  large  t  or  very  small  time  constants,  the  OC  approaches  classical 
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Appendix  B 

NUMERICAL  DATA  FOR  OPTIMUM  CONTROLLER  GA:NS 

The  numerical  values  for  the  functions  N,  G:i  and  G4  were  computed  via 
digital  program  for  a  wide  range  of  values  of  w,  and  u2.  The  data  are  presented 

in  this  appendix  in  the  form  of  plots  of  the  subject  functions  versus  t  ^Figures 

go  H 

B-l  through  13-19).  The  time  span  considered  is  10  seconds,  which  should 
cover  most  practical  applications.  The  assymptotie  nature  of  the  functions 
allows  easy  extrapolation  from  the  plots  if  somewhat  ionger  time  spans  are 
deemed  desirable.  In  some  instances,  particularly  where  u,=  u2,  the 
computation  of  N  involves  differences  of  large,  almost  equal  numbers;  for  this 
reason,  double  precision  arithmetic  was  found  necessary. 
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